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1 . Trianyulation 

2. Trilateration 

3. Traverse 

4. Intersection 

5. Resection 

The process of making a proper nautical chart consists, 
first of all, in setting up a framework of marks on the 
ground- Before 1950 the main framework of a geodetic survey 
almost always consisted of triangulation, which was replaced 
by traverse if the topography made triangula tion impracti- 
cable. During the last decade, the introduction of elec- 
tronic distance measuring (EDM) equipment has made both 
trilateration and traverse economical, and an acceptable 
substitute for triangulation. In fact, it appears probable 
that these new methods will replace triangulation as the 
main framework for new geodetic surveys [Ref. 1]. It is 
evident that triangulation and traverse are the main methods 
used for establishing control- There seems to be no 



lA position of a point on the surface of the Earth 
expressed in terms of geodetic latitude, geodetic longitude, 
and geodetic height. A geodetic position implies an adopted 
geodetic datum. 
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agreement among the various agencies as to which of these 
two methods is mostly used. The Q.S. tiational Ocean Service 
(NOS) does the majority of its horizontal control surveys 
for hydrography with traverse (about 90 %) [Eef. 2]. The 
main factor for the selection of one or the other method 
depends on the geographical configuration of the survey 
project area and the availability of good EDM equipment. 
Traversing is a convenient, rapid method for establishing 
horizontal control. It is particularly useful in densely 
built areas, when the coastline tends to be even, along a 
railroad track, and in heavily forested areas where lengths 
of sight are short so that neither triangulation nor 
trilateration is suitable. 

The objectives of this thesis are to show methods of 
traverse station position computation and a comparison of 
methods of least squares adjustment. Closed traverse 
station positions were computed and adjusted using the 
Approximate Method and by the Least Squares Method in the 
Oniversal Transverse Mercator (UTM) grid coordinates. The 
computer programs were written to calculate the adjusted 
traverse station positions. Test data included those data 
obtained during the Geodetic Survey Field Experience course 
at the Naval Postgraduate School (NPS) in October 1983. 

The results of comparative computations are shown to 
more significant figures, in this thesis, than are normally 
considered desirable in production work. The same observed 
data are used with each of three computational methods. It 
is important to recognize that what is being compared is not 
observational precision but computational precision. Hence, 
it is considered necessary for a rigorous comparison of 
computational precision, including round off error, to show 
results to several more decimal places than is justifiable 
based on observational precision alone. 
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TRAVERSE 



A. GENERAL 

Traverse is a method of surveying in which a series of 
straight lines connect successive established points along 
the route of a survey. An angular measurement is taken 
using a theodolite at each point where the traverse changes 
direction. Distances along the line between successive 
traverse points are determined by EDcl equipment. The points 
defining the ends of the traverse lines are called turning 
points, traverse points, or traverse stations. Each 
straight section of a traverse is called a leg or a traverse 
line . 

A closed traverse originates at a point of known posi- 
tion and is closed or. another point of known horizontal 
position. Traverse I-1-2-3-F originates at point I with a 
backsight along line lA of known azimuth and closes on point 
F, with a foresight along line FB also of known azimuth 
(Figure 2.1). This type of traverse is preferable to all 




Figure 2. 1 
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Figure 2.2 Closed-loop Traverse. 



B. ANGLE AND DIRECTION MEASOREMENT 

Angles and directions may be defined by means of b 
ings, azimuths, deflection angles, angles to the right 
interior angles. These quantities are said to be obse 



^Systematic errors may be caused by faulty instrum 
or factors such as temperature or humility changes whi 
affect the performance of measuring instruments. 
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when obtained directly in the field and calculated when 
obtained indirectly by computation. 

A theodolite is an instrument designed to observe hori- 
zontal directions and measure vertical angles. It consists 
of a telescope mounted to rotate vertically on a horizontal 
axis supported by a pair of vertical standards attached to a 
revolvable circular plate containing a graduated circle for 
observing horizontal directions. Another graduated arc is 
attached to one standard so that vertical angles can be 
measured . 

Repeating and direction theodolites have features that 
are common to both types of instruments. Repeating theodol- 
ites are read directly to 20" or 01’ and by estimation to 
one-tenth the corresponding direct reading. Direction theo- 
dolites are usually read directly to 01" and can be esti- 
mated to tenths of seconds [Ref. 3, p- 215]. In general, 
direction theodolites are more precise than are repeating 
theod elites. 

The direction theodolite observes directions only and 
angles are computed by subtracting one direction from 
another. Assume that the horizontal angle AI3 (Figure 2.3) 
is to be measured with a direction theodolite. The theo- 
dolite is set over point I, leveled and centered, and a 
sight is taken on point A. The horizontal circle is then 
viewed through the optical- viewing system and the circle 
reading is observed and recorded. Assume the reading is 
45° 02' 40". The telescope is then sighted on point B. The 

horizontal circle is then viewed through the optical-viewing 
system and the circle reading is observed and recordel. 
Assume the reading is 124° 11' 59". These two obs'ervations 
constitute directions which have a common reference direc- 
tion that is completely arbitrary. The clockwise horizontal 
angle is 

i = (1240 in 59«i) - (450 02 ’ 40") = 79° 09' 19" 
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Figure 2.3 Horizontal Angle. 



C. DISTANCE HEASOREHENTS 



There are several methods of determining lis 
choice of which depends on the accuracy required 
cost. For example, tacheometry, taping, and ED.1 
can be used. The general availability of ED!1 eg 
practically eliminated the use of taping for the 
of traverse lengths. Accuracies are comparable, 
superior, to those obtained with invar tapes. 

Distances measured using EDM equipment are s 
errors arising from the instrumental components, 
of the equipment, inaccuracies in the meteorolog 
elevation discrepancies , and centering of the in 
reflectors. The reduction of measured distances 
converting the slope distance to a horizontal di 
converting the chord distance to an equivalent a 
and reducing the arc distance to the ellipsoid, 
tion of slope distance to horizontal distance is 
to compensate for the difference in elevation of 
points of the measured line. The horizontal dis 
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TABLE II 

Coor3inates of Known Stations and Aziauths 



T 



Station 


UTM grid 
Grid Northing 
( m. ) 


coordinates 

Grid Easting 
( m. ) 


Moss 2 


4, 


072,555. 85206 


608,279.04404 


Holm 


4, 


079,258. 31754 


612, 238.85256 


Grid 


azimuth clockwise 
from North 


Azimuth 


From Moss 2 


to Pipher 


s tations 


0 » n 

100 16 23.778 


From Holm 


to Moran 


stations 


136 33 26.334 



determination. The distances were observed in the field, 
corrected by temperature and pressure for propagation error. 
Unfortunately, meteorological data are generally acquired 
only at the end points of a measured line. Using the mean 
of these meteorological values only approximates the actual 
conditions of the entire measured line and does not 
completely correct for errors in the propagation velocity of 
electromagnetic radiation. By applying the elevation, sea 
level, and scale factor corrections they were reduced to UTH 
grid distances [Ref. 4, pp. 124- 125 ]. The UTM grid 
distances and standard deviations of the distances were 
determined in meters (Table III) . 

The angles were observed by using a Wild T-2 theodolite. 
To ensure the correctness of the beginning and ending 
azimuths, a check azimuth to a second station of known posi- 
tion was observed. The angles were observed at stations 
Moss 2, Mossback, Dune Temp, and Holm (Table IV). All 
observations were made by NPS students and conform to 
specification for a third-order class I traverse. 



TABLE III 

Grid Distances and Standard Deviations 



I 1 

Standard 

Grid distances between Distances deviation 

stations ( m. ) ( m. ) 

Hoss 2 and !1ossback 1/424.004 0.00 1 

Mossback and Dune Temp 365.744 0.001 

Dune Temp and Holm 6,476.271 0.003 

TABLE I? 

The Observed Angles and Standard Deviations 

Backsight Center Foresight Observed Standard 

station station station angles deviation 

0 • I* 

Pipher Hoss 2 Mossback 246 05 43.200 01.984 

Hoss 2 Mossback Dune Temp 222 51 08.600 01.405 

Mossback Dune Temp Holm 190 15 02.600 01.203 

Dune Temp Holm Moran 277 05 17.000 01.614 
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III. TRA7ERSE COMPDTATIONS ASD ADJOSTMSNTS 



A. INITIAL TRAVERSE COMPOTATIONS 

A traverse (Figure 3. 1) originates at station 1 (known 

position) and terminates at station 4 (also known position) 

(Table V). To compute the forward azimuth (Table VI) of an 

unknown leg, the angle is added to the back azimuth of the 

previous leg (Equation 3. 1) . 

i 

AZj = Az 2 “ (1-1)180° (3.1) 

j : 1 

Where: i = the number of legs, 

AZj = the forward azimuth of an unknown leg, and 
AZp. = the fixed initial azimuth. 

For the computation of OTM grid coordinates, let AEj 
and ANj be designated as the departure and latitude for leg 
i (Figure 3.1) so that the general formulas to compute the 
departure and latitude are 



It 

w 

0 


d; 


sin AZj 


(3.2) 


0 

H 


dj 


cos Az. 

1 


(3.3) 


where 


i 


— 1,2,3,.. .n. 





n = the number of the observed angles. 

The algebraic signs of the departure and latitude for a 
traverse leg depend on the signs of the sine and cosine of 
the azimuth of that leg. The algebraic sign of the depar- 
ture and latitude is determined by the following rules: 

1. DTM grid azimuths are refered to grid north, 

2. For azimuths between 0° and 180°, the departure is 
plus; for all other azimuths, the departure is minus. 

3. For azimuths between 90° and 270°, the latitude is 
minus; for all other azimuths, the latitude is plus. 
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Figure 3. 1 A Closed Traverse. 



The calculators and computers with internal routines for 
trigonometric functions yield the proper sign automatically 
given the azimuth. The major portion of traverse computa- 
tions consists of calculating departures and latitudes for 
successive legs and cumulating the values to determine the 
coordinates for consecutive traverse stations. 
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TABLE V 

Data for Initial Traverse Computations 



Leg 


Distance 




Angle 


( m. ) 


^1 


0 


f n 


II 

CM 

1 


1,424, 004 


= 246 


05 43.200 


2-3 ^2 = 


365. 744 


C<2 


= 222 


51 08.600 


3-4 ^3 = 


6,476.271 




= 190 


15 02. 600 


The initial OTM 


grid azimuth 


AZp 


= 100 


16 23. 778 


QTM 

Station 1 


grid Northing 




= 4,072 


,555.85206 m 


DTK 


grid Easting 


El 


= 608 


,279.04404 m 



The coordinates of station j (j = i 1) are Ej and 
Nj , the coordinates of station j are 

Ej = E. + AE, (3.4) 

N. = Hj + AN| (3.5) 

wherei=1 , 2 , 3 , . . .n. 

The values of departures, latitudes, and coordinates are 
shown in Table 7II and were computed from the data in Tables 
V and VI by application of Equations (3.2) through (3,5). 

B. COMPOTATION OF DISCREPANCIES 

A closed traverse at Moss Landing (Figure 3. 1) origi- 
nated at station 1 (Moss 2) and closed at station 4 (Holm). 
The closing errors for a traverse are caused by observation 
errors in the observed angles and the measured distances. 

The closing errors may be computed by applying Equations 
(3.1) through (3.6). 
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TABLE VI 

The Azimuth Calculation 



— 
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F or ward 
Azimuth 


Back 

Azimuth 


Observed 

Angle 


0 




It 






AZp 




100 


16 


23.778 








c<i 


246 


05 


43. 200 




Az^ 






346 


22 


06. 978 
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180 


00 


00. 000 








^2 


166 


22 


06. 978 








222 


51 


08.600 




Az2 






389 


13 


1 5. 578 








- 


180 


00 


00.000 






AZ3.2 




209 


13 


1 5. 578 








190 


15 


02. 600 




A 






399 


28 


18. 1 78 








- 


180 


00 


00. 000 










219 


28 


18.178 














J 



The angular error of closure may be computed by applying 
Equation (3.1) becomes 

n 

Az + 2ci.- (n-1) 1800 - Az = (3.6) 

i 

where Az^^ is the fixed closing azimuth and is the 
angular error of closure. Az^ is the grid azimuth from 
stations Holm to Horan (AZj^ = 136 ® 33' 26. 334") (Table II). 
The computed closing azimuth is equal to Az^_^ plus the 
observed angle at station Holm = ( 219 ® 28' 18.178") + 

(2770 05' 17") = 496 0 33’ 35. 178" = 136° 33 * 35. 178". From 

Equation (3.6) , the angular error of closure is 
(1360 33 ' 35. 178") - (1 360 33' 26.334") = + 08. 844". 
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TABLE 7II 

Initial Traverse Computations 



i — 




Values 


in meters 






Values in meters 




Ni 


4,072, 


555.85206 




Ei 


608,279. 04404 


An, 






383.89267 


AE, 




- 335. 60166 


'^2 


4,073, 


939.74473 


E 2 


607,943.44238 


AN 2 






31 9.20061 


AE 2 




178.54871 




4,074, 


258.94534 


E 3 


608, 121.99109 


AN 3 




4, 


999.28234 


AE 3 




4,116 . 94755 


'’4 


4,079, 


258.22818 


E 4 


612,238. 93864 
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.5) are applied to 


obtain; 
















n -1 












?2 = 


2 ae, 

i : 1 


- ( Et - 


El ) - 




(3. 7) 






n -1 












«3 = 




- ( - 


N, ) 




(3.8) 


or 




i r 1 












^2 = 


En - 








(3.7a) 




"3 = 


Nn - 








(3.8a) 



Where: n = the number of observed angles, 

W '2 , W 3 = precalculated discrepancies, 

= the fixed initial coordinates, 

E-j. , = the fixed closing coordinates, and 

, Np = the computed closing coordinates. 

The result of the errors in closure in position of a 
traverse is 
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(3.9) 



Ad = ( W I + H I ) 1 / 2 

The computed closing coordinates of this traverse are 
= 612,238.93 864 m and = 4,079,258.228 18 m. The fixed 
closing coordinates are Ej = 612,238.85256 m and = 
4,079,258.31754 m. The precalculated discrepancies can be 
computed by using Equations (3.7a) and (3.8a). The results 
are W 2 = 0.08608 m and = -0.08936 m. 3y applying 
Equation (3.9), Ad is 0.12406 m. From Table III, the total 
measured distance is 8,266.019 m. The ratio of the distance 
error of closure to the total distance is 0.12408/3,266.019 
or 1 part in 66,618. The ratio is an indication of the 
goodness of this traverse. 

C. ADJOSTHENT OF A TRAVERSE BY AN APPROXIMATE METHOD 

The angular error of closure of traverse may be distrib- 
uted equally among the observed angles. The angular error 
of closure of this traverse is +8. 344'*, which corresponds to 
a correction for each angle of -2.211”. The observed angles 
were corrected by this value. The adjusted azimuth was then 
computed. The departures, latitudes, and coordinates were 
recomputed by using the adjusted azimuth (Table VIII). The 
closure corrections in E and N coordinates are +0.09635 m 
and -0.04326 m, respectively. The resultant closure is 
0.10562 m. The ratio of the distance error of closure to 
the total distance is 0.10562/8,266.019 or 1 part in 78,262. 

The adjustments of a closed traverse by the Approximate 
Method are completed by the compass rule which proportions 
the errors in E and N coordinates according to the distance 
of the course [Ref. 4, p. 354 ], The corrections are applied 
to the departures and latitudes prior to computation of 
coordinates. 



\ 
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TABLE VIII 



Initial Traverse Computations (Approximate Method with 

Adjusted Aziamtn) 





Observed angle Correction 


Adjusted 


azimuth 








0 


f 


II 






A Z|_ 


100 


16 


23.778 




2460 


05’ 43.2" - 2.211" 


246 


05 


40.989 






Azi 


346 


22 


04. 767 








- 180 


00 


00.000 






A z 2-1 


166 


22 


04. 767 


“^2 


2220 


51' 08.6" - 2.211" 


222 


51 


06.389 






A Z2 


389 


13 


11.156 








- 180 


00 


00. 000 






AZ3_2 


209 


1 3 


1 1. 156 




1900 


15' 02.6" - 2.211" 


190 


15 


00. 389 






A Z^ 


399 


28 


1 1 . 5 45 








- 180 


00 


00. 000 






A z, ^ 


219 


28 


1 1 .545 






4-3 










2 77 0 


05 ' 17.0" - 2. 21 1 " 


277 


05 


14. 789 






AZ 4 -L 


136 


33 


26.334 






Values in meters 


Values 


in 


meters 




Ni 


4,072,555.85206 S, 


608, 


279. 


04404 


ANi 




1,383.88907 AE^ 


- 


335. 


61649 




^2 


4,073,939.74113 S2 


607, 


943. 


42755 


An, 




319.20444 AE. 




178. 


54137 






4,074, 258-. 945S7 E3 


608 , 


121. 


96942 


A N. 




4,999.41523 AE. 


4, 


116. 


78679 




«4 


4,079,258.36080 E^ 


612, 


238. 


75621 






4,079,258.31754 E^ 


612, 


238. 


85256 




1 

dN 


0.04326 dE 




- 0. 


V 

09635 
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(3.10) 



Where : 



5 E. = ( dE / D ) . d [ 

5N| = ( dN / D ) . d. (3.11) 

5Ej = correction to 
5 Nj = correction to AN] , 

dE = total closure correction in the E coordinate, 

dN = total closure correction in the N coordinate, 

and 

D = total distance. 



The adjusted S and N coordinates (Table IX) were 
different from the fixed E and N coordinates (Table II) by 
0.00001 m due to round off error. The calculations and 
adjustments were illustrated in sections A, B, and C by 
using the hand calculator (TI-59) and rounded off at 5 
decimal places. The computer program was written to calcu- 
late the adjusted traverse station positions by the 
Approximate Method by using 16 decimal places (Appendix A). 
The approximate traverse adjustment is based on an assumed 
condition that the angular precision equals the precision in 
linear distance. 



D. ADJUSTMENT OF A TRAVERSE BY LEAST SQUARES METHODS 

The method of least squares provides a rigorous 
adjustment and best estimates for positions of all traverse 
stations. The Least Squares Method is used to 
simultaneously eliminate closing errors in azimuths and 
coordinates of traverses, 

"I • The Princi pi e of Least Sguares 

The fundamental condition of the least squares tech- 
nique in surveying requires that the sum of the squares of 
the residuals be minimized. A residual is defined as the 
difference between the true and observed values. In making 
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TABLE IX 

Adjusted Coordinates (Approximate Method) 



Departure Correction 
AE, -335.61649 + 0.01 659 
AE 2 178.54187 + 0.00426 
AE 3 4,1 16.78679 + 0.07549 

Latitude Correction 



AN, 


1, 383. 88907 


- 0.00745 


AN2 


3 19 . 20444 


- 0.00191 


AN3 


4,999. 41523 


- 0.03389 



Grid Values in meters 
608,279.04404 
- 335.59990 
E 2 607,943.44414 

178.54613 
E 3 608,121.99027 

4, 1 16. 86228 
E 4 612, 238. 85255 

Grid Values in meters 
4,072,555.85206 
1, 383. 38162 
N 2 4,073,939.73368 

319.20253 

N 3 4,074,258.93621 

4, 999. 38134 
N 4,079,258.31755 

4 . 



physical measurements, the true values can never be deter- 
mined. The least squares principle establishes a criterion 
for obtaining the best estimates of the true values. 

If the best estimates of the true values are stated 
by X. and observed values by x. , the residuals are expressed 
as 

V, = X, -X,. 

The fundamental condition of least squares, for 
uncorrelated observations with equal precision are expressed 
as 
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n 

2 ^ i ) ^ = (^1 ) ^ ^ (^ 3 ) ^ • * * ■•■ ^ " Minimum 

i : 1 

or in matrix form 

v''’v = Minimum (3. 12) 



In general, the observed values are of unequal 
precision. The observed value of high precision has a small 
variance. Conversely, a low precision of the observed value 
has a large variance. Since the value of the variance goes 
in opposite direction to that of the precision, the observa- 
tion is assigned a value called weight corresponding to a 
quality that is inversely proportional to the observation’s 
variance. 

For uncorrelated measurements x, , x-, x^ , . . . 

^222 2 ^ 

Xp , with variances 0 ^, ... , respectively, the 

weights of these uncorrelated measurements are, 

p = ^ 0 / P = CL, p = 6 . 

11 ^22 2 33 3 



2 2 

Pnn= (3.13) 

n n 

2 

where O’, is the proportionally constant of an observation 
of unit weight [Ref. 7, p. 67]. These weights may be 
collected into a corresponding diagonal P matrix, called the 
weight matrix: 



P = 




0 

0 

0 









0 



0 



0 



p 



n n 
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The weight matrix consists of weighted o 
of a traverse with mixed kinds of measurements, 
distances and angles- The variance of a measure 
is expressed in meters and an angle in radians. 
uncorrelated observations of unequal precision, 
mental condition of least squares is expressed a 




p + p v 2 + p 72 

^11 1 ^22 2 ^33 3 

Minimum 



or in matrix form 

V^PV = Minimum 



f 



valu 

of t 

rela 

rela 

expa 

‘or de 

solv 

the 

mini 



In the problems involving the adjustment 
es, all observed quantities are expressed by 
he quantities to be determined. In simple c 
tions are linear, but when this is not the c 
tions must’ be converted into the linear form 
Tiding them into Taylor’s series. The terms 
r are neglected, so as to obtain linear rela 
ing the resulting linear equations, then ita 
effect of the neglected higher order terms a 
mixed. 



2. least Squares Ad jus t me nt of Indirect Obs 

The observed values are related to the d 
unknown values through formulas or functions whi 
called observation equations. Cne observation a 
written for each measurement. To solve for the 
of each unknown parameter, at least one redundan 
tion equation must be written. That is, the num 
tions must be greater than the number of unknown 
linear observation equations can be written in t 
form as follows: 



3 1 



bser vat ions 
They are 
d distance 
For 

the funda- 
s 

p V 2 
n n n 



(3.14) 

of observed 
f unct ion s 
ases these 
ase, the 
by 

of higher 
tions, 

rating until 
re 

er v a ti ons 

esired 
ch are 
quation is 
best value 
t observa- 
ber observa- 
s. The 
he general 



10 


^11^1 


^ a^2^^2 ^ * * * 


* 


= 


Gi + v^ 


20 "^ 


^21^1 


+ a22^2 ■*■••• 


* ^2m 


= 


G 2 + V 2 






^ ^n2^2 ^ 



+ a X = G + 7 
n m m n n 



where n is the number of observations; m is the number of 
the unknowns; a^Q, ^20' ' * * ^no constants; a^^ , ^12' 

. . , ^nm coefficients of the unknowns , x^t • • • 

x^ ; and » ^2 ' * * * residuals. 

Because the observations Gj (i = 1, 2, . . . n) are 

not free from random errors, each Gj must be corrected by a 
residual value, Vj . let b- = Gj - S|q- Thus 



1 1 


^1 


^12^2 ^ * * 


. ^ 


^lm 


- 


bi = v^ 


2 1 


^1 


+ ^22^2 • 


• ^ 


^2 m 


- 


b2 = V 2 



^nl^l 



^n2^2 



a X 
n m m 



b 

n 



V 

n 



or in the matrix form 



V = AX - B 



(3.15) 



This equation is called the observation equation or 
the observation equation matrix. 

For uncorrelated observations of unequal precision, 
substitute the value for the 7 matrix from the observation 
Equation (3.15) into Equation (3.14) 
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V^PV = (AX - B)"'’p(AX - B) 

= [ (AX)'*' - B^]P(AX - B) 

= (X’’’ a"^ - b’’’) P (AX - B) 

= (x'’’a''’p - b''’p) (AX - B) 

= x''’at PAX - x'^a'’’pb - b''’pax + b'^’pb 

= x''’a^pax - x''’a''’pb - x'''a''’pb + b''’pb 

= x''’a"pAX - 2x'^a‘'‘pB + b''’pB (3.16) 

[from matrix algebra, (AX)"'" = x''' a"*" and b"''paX = x'*'a'^PB] 

The minimum of this function can be found by taking 
the partial derivatives of the function with respect to each 
unknown variable (i.e., with respect to the x^ , X 2 » • • • 
x^ ) must equal zero. Hence, 

^(vTpyj ^ 2 a''’ PAX - 2a"^PB = 0 (3.17) 

dx 

Dividing Equation (3.17) by 2, the following result 
is obtained: 

a’'’pax = a''’pb 

This represents the normal equations, and by multi- 
plying this equation by (A^PA)~i, the solution is obtained: 

X = (a''’pA) -ia''’pb (3.18) 

For uncorrelated observations with equal precision, 
the weight matrix is an identity matrix, and equation (3.18) 
becomes 

X = (A^A)-ia‘''b (3.19) 

This equation can be derived similarly to the 
unequal precision case. Equations (3-18) and (3.19) are the 
basic least squares matrix equations. 
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If the relations are nonlinear, the relations must 
be converted into the linear form by using Taylor's series 
expansion [Eef. 3, p. 919], Let F = f (x^ , • • • x ) 

be the general observation equation that is a nonlinear 
function. The Taylor's series expansion is 

F = f (2* , , - . X° ) + 5 2, ♦ ^ 5 2j* - • * Ax 5 ^rn 

3x- ax-2 

+ Higher order terms 



where x° , x° 2 / • • - are the approximate values of the 

variables at which the function is evaluated. For the 
traverse problems, an approximate value can be the 
precalculated value, and 5x^, 8 x 2 , . . . S x^^ can be the 
corrections. The higher order terms in the series are 
neglected and only the zero and first order terms are 
maintained, the approximate value must be improved by 
successive iterations until the effect of the neglected 
higher order terms is minimized. After linearization, the 
observation equations become: 



, = f, (x“, X 2 , . . ^x, + ^,5x2 + . . +^i5x^-G, 



3x 



m 



= f2(x“/ x°2/ . - X°^) + 






m 






2' 




IfnSxi ^n^X2 + 
3x°2 



or in the matrix form 
V = AX - B 






■ m 



(3.20) 
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wher e 











< ^ 


'dx\ 


'bty . . . 

'3x°2 








U 2 


3f, . . . 

3x*2 


. , "3 f2 

'dx 
^ m 




^2 


- 


m 




V = 


• 


• 








m 


• 


- 


• 




m 


Un 




. . 3fn 






"dx' 
^ ' 


d>^2 


^x° 
m J 

















Gj “ fj (x.| , x^/ .... ) 




5 




^2 ■" £ 2 ^^!' ^2* ’ * * * 




6., 


B = 


« • • • • 


X = 


• 




« « « • • 

^ ~ (^1 » ^ 2 ' * * * * ^ ^ 




6 X 

mj 



The remainder of the least squares procedure is the 

same as indicated by Equations (3.16) to (3.19). 

There are two types of observation equations for the 

adjustment of the traverse. They are the angle and distance 

observation equations. Both of them are nonlinear, and they 

must be linearized using the Taylor's series. 

To derive the observation equation for the angle c(. 

(Figure 3.1), first write the two azimuth observation 

equations for the azimuth Az-j^ and AZjj , where Az.^^ is the 

back azimuth and Az.. is the forward azimuth. 

ij 

Az,^ = arc tan = 

^j-^i 

Az.. = arc tan m. ■ _ m — 

IJ 
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Then, the angle observation equation for the angle c( • is 



V. = Az. . - Az., - Ol. 

I I j I k I 



= arc tan 



£i-^i 

Nj - N| 



^ Ek-Ei 
- arc tan 



N^-Nj 



- d: 



(3.21) 



Equation (3.21) is nonlinear in the parajieters. 
Write this equation in function form: 

■ 7 . = £, (Z;, H;, Ej, N,, E^, S, ) - d, 



Thus, the 


li nearized 


form of 


the angle 


observation 


equation is 










V, = f; (E] , N 


0 0 0 

i ' ^}' ' 


0 0 


+ B f . 5E; 


+ B f : 5 N • 




k ^ 


Be] 


Bn] 


+ B f. 5e- 


+ B f, 8N; 


+ B f ; 5 E , 


+ B f ; 5 - 




' 






Bn“, 

k 


The observation equation for 


the distan 


ce dj between 



two points i and j is 

V, = C (E; - Ej) 2 + (Nf - Nj ) 2]i/2 - d 

or in the functional form 



V; = f i (S| , N, , Ej, Nj ) - d 



Thus, the linearized form of the distance 
observation equation is 

V; = f;(E;, Nj, E° N^- ) + f 5 E , ^ ^;5N; 



Be 



Bn: 



+ 9f. 5e- "Sf-SN- -d- 

■ ■ ■ I I — — = I I I 

Be] Bn] ^ 



For the adjustment of the traverse which was 
conducted at tioss Landing with the technique of adjustment 
of indirect observations, four angle and three distance 
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observation equations need to be written. The seven equa- 
tions would include as unknown paracieters the four coordi- 
nates of stations 2 (Mossback) and 3 (Dune Temp) . Four 
normal equations are formed and solved for corrections to 
the approximate values (precalculated values) of the parame- 
ters. The corrections are added to the approximations to 
update their values, and the solution is repeated until the 
last set of corrections is insignificantly small. 

The linearized equations of this problems are 



v^ = a 



V2 = a 



^7 = ^ 



'■1 1 ^^2 


♦ a„S», 


+ a^3 6E3 ♦ a, 4 5 n 


i 5 e 


<• a^^SN^ 


a^^SE., + a^^^N 


21 2 

• 


22 2 

m 

• 


2 3 3 2 4 

« • 

• « 


« 


• 

+ a S N 

72 2 


• • 

* ^73 5E3 ♦ 6s 



- b- 



The zero order terms of the angle and distance 
functions are 

E o - E 



fi = N=2 -Ni 

f2 = arc tan ^ — — 



- Az, 



- arc tan 



f3 = arc tan ^ ^ 



N|-N^ --- Ni-n°2 



- arc tan 



El - E2 

Nl - N°2 
^ 2-^3 

- N »3 



E^- E. 

^4 = n|-N, 

3 4 

f^ = [ (E^ - E2 ) 2 ♦ (N^ - N2 ) 2 ]t /2 



fg = [ (E 2 - E 3 )2 + (N 2 - N °3 ) 2] 



2 11/2 



f ^ = [ (E3 - E^ ) 2 + (N! - N4 ) 2 ] 



1 /2 



By using the data of the known coordinates (Table II) and 
the approximate values (Table VII) 



a 



11 



= 

3 e »2 



N2 - Ni 

- -I- (E2- El )^_ 



0. 68246 X 10-3 
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0. 16550 X 10-3 



^12 = 



M 3 






E? - S- 



l(N 2-N^)2 +- (E*- )^j 



= 0 



^’4 =1^1 



= 0 



‘21 



= lf = - 


N3 - N°2 




Ni - N 2 j 


'Se-2 


£N3-N2 +- (E°3- £“2 )^. 




iN^-N2)2+ (Et-E°2)2. 



= - 3. 06868 X 10-3 



‘22 



= 3 f 2 = 


r~ 00 ^ 

E 3 ~ E 2 


- 


El - E 2 


■^N °2 


_CN 3 ~ N *2 ) ^ + (E 3 - E° 2 )^_ 




W^~ N2)2 +. (E^- E° 2 )-_ 



= 1. 16926 X 10-3 

O o 

M-:; - N? 



^23 =^2 



3e° 



^24 ~ ^ 2 



^ N 



i^3~ N2)2-^ (E;-E“)' 



E3 - Ea . 



3 |(N-3-Np2 ^-<e-3-Ej)-^J 



^31 

■^Eo 



^3 2 -3-i3 



N2 - Ni 



9n° 



[(N2-N3)2^(E“2-Ey 

E2 - 



1_(N2-N3)2 +(E2-E°3)2 



= 2.38621 X 10-3 



= - 1.33476 X 10-3 



= 2;38621 x 10~3 



= - 1.33476 I 10-3 



‘3 3 



= 3 f 3 = - 


r N4 - n; 


> 


r 0 0 ^ 

N2 - N3 


^^3 


[(N^NjJSmE^- E-3 )2^ 




N°3)2 ^ (E»_ £»^)2^ 



= - 2.50541 X 10-3 



3 4 :^3 

^N<’3 





E4 - E3 


- 


r 0 —0 ^ 

E2 - E3 




iN 4 -N 3)2 t-(E^- E^)2J 




jN^- n!^)2 ^(E2 -E3)2 



= 1.43291 X 10-3 
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0 



a , = 9 f ^ 

4 1 ■::: — ^ 

dE°2 



^2 = M. ' » 

3n'’2 



a.o - ^ 

4 3 ^57:4 



a = S 

4 4 

9n°2 



N^ - N 4 



\ 1 ^ 3 ~ ^ 4 )^+ (E“-E^) 2 J 

~ E 4 



^N 3 ~ N^)2 ^ (E 3 - E4)2j 



= 0. 1 1920 X 10-3 



= - 0.09816 X 10-3 



a = 9 f 

51 -^5 

9e^2 



a = 2 =-^5 
SNf 



El - S2 



_[(E,~ E^)^ 4- (N^ — N2)2]^^2 

Ni - N2 



f(E,- E°2)^4 -(Nt - N“2)2]^'2 



= - 235,67466 x 



= 971.83201 X 10 



a-3 =^£5 

9e»3 



= 0 



= 0 

5 4 D 

3n’3 



0£% 



ag 2 =-^6 
9Nf 



^63 "—6 
9e’3 



a = ^fg 

"JSf 



■* 


0 0 

E? - E3 


“ 


.[(E’-E’^)^ -|-(N2-N3 )2 ] 


1 /2 


r 


0 0 
N 2 - N 3 


- 


_[(E2-E;)2 +-(n“2-N“3)2] 


|1 /2 


- 


B 2 - E 3 


- 




JjE2-E3)^ -h(N°-N^)2] 


1 /2 


- 


N °2 - N 3 


-1 




1 (E 2 - E3)2 +(N2~ N°3)2; 


jl /2 



= - 488, 1 7948 x 



= - 872.74326 x 



= 488.17948 x 10 



= 872.74326 x 10 



a,, =11,7 = 0 

3E2 



a^2 =-ify = 0 

^ 9n2 



10-3 

—3 

10-3 
10 -3 

-3 

-3 
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635.68254 x 10 "3 



73 



” ^"7 " 



E3 - B, 



J(E^- E^)2^(N°3-N^) 



2 11/2 



■74 



= 3f, = 



- N. 



I(E3-E^) +(N3-N^) 



2 | 1/2 



= - 771. 95059 X 10-3 



or in the m 



A = 10-3 



rix form 

0. 68246 
-3. 06868 
2.38621 
0 . 00000 
-235. 67466 
-448. 17948 
0 . 00000 



0.16550 
1. 16926 
-1.33476 
0.00000 
971.83201 
-872.74326 
0. 00000 



0.00000 
2.38621 
-2.50541 
0. 11920 
0. 00000 
488.17948 
-635. 68254 



0 - 00000 
-1.33476 
1. 43291 
-0. 09816 
0 . 00000 
872.74326 
-771. 95059 



The computation of B matrix, 
from Tables II, III, IV, and VII 



by using the data 





f^ = 0.00000 radians 

f 2 = 0.00000 radians 

f 3 = 0.00002 radians 

f^ = 0.00002 radians 

f^ = 0.00000 meters 

fg = 0.00000 meters 

fy = -0.01426 meters 



or in the matrix form 

= [ 0 0 0.00002 0.00002 0 0 -0.01426 ] 

The weight matrix can be computed by using the vari- 
ance of the observed angles and distances from Tables III 
and IV by applying to Equation (3.13). For uncorrelated 
observations of unequal precision, the weight matrix is 
defined as the diagonal ? matrix. 
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P 11 


= 


1 


/ 


sin 2 (1.984") 


= 


10,808,537,426. 79957 


?2 2 

P33 




1 


/ 


sin 2 (1.405") 


= 


21,552,498,219.02474 


= 


1 


/ 


sin2 (1.203") 


= 


29,398,082,997.43487 


P 


= 


1 


/ 


sin 2 (1.614") 




16,332,144, 194.08730 
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^5 5 

p 


= 


1 


/ 


(0. 001) 2 


— 


1,000,000.00000 




1 


/ 


(0. 001) 2 


= 


1,000,000. 00000 


•66 


^77 




1 


/ 


(0. 003) 2 


- 


111,111.11111 



For uncorrelated observations of unequal precision. 



the values of A, B, and P matrices are applied to Equation 



( 3 . 18 ). 



The correction vector is found to be 



X = 



5e2 

5n^ 

5e3 



= 0.01284 
= 0.00336 
= 0.01079 
= 0.00495 



The approximation values were improved by adding the 
corrections to the first approximation values. The improved 
approximation values after the first iteration are 



E2 = 607,943.44238 + 0.01284 = 607,943.45522 m. 
N2 = 4,073,939.74473 + 0.00336 = 4,073,939.74809 a. 
E3 = 608,121.99110 + 0.01079 = 608,122.00189 m. 
N3 = 4,074,258.94534 + 0.00495 = 4,074,258.95029 m. 



Using these values, the solution is iterated. After the 
second iteration, the correction vector is zero to six 
decimal places, so the improved approximation values are 
the final estimates of the coordinates. 

For uncorrelated observations of equal precision, 
the correction vector is found by solving Equation (3.19) . 
The correction vector after the first iteration is 



5 E2 


= 0.01756 




= 0.00426 


5E3 


= 0.01660 


5 H3 


= 0.00479 



The imprcved approximation values after the first iteration 



are 



E 2 = 607,943.44238 + 0.01756 = 607,943.45994 m. 

N 2 = 4,073,939.74473 + 0.00426 = 4,073,939.74899 a. 

E 3 = 608,121.99110 + 0.01660 = 608,122.00770 m. 

N 3 = 4,074,258.94534 + 0.00479 = 4,074,258.95013 m. 

The solution is iterated by using these adjusted 
values. The correction vector is zero to six decimal places 
after the second iteration. These values are the final 
estimates of the coordinates. 

The standard deviation cf an observation vhich has 
unit weight can be found by the following equations, for 
uncorrelated observations of unequal precision. 



( 5 „ = 



P7 

n - m 



X / 2 



(3.22) 



for uncorrelated observations of equal precision. 



d. = 






1 /a 



(3.23) 



n -m 

where n is the number of observation equations and m is the 
number of unknowns [ Eef . 7, p. 249]. After calculating the 
best estimate values of the unknowns, the V matrix can be 
computed from Equation (3.20). 

The standard deviations of the best estimate values 
for the unknowns are then given by the following equations: 




X / 2 



(3.24) 



where C- is the standard deviation of the ith adjusted 
quantity. The quantity in the ith row of the X matrix, S- 
is an element of the (a"*^ PA) -1 matrix for uncorcelated 
observations of unequal precision. For uncorrelated 
observations of equal precision, Sj] is an element of the 
(a"’’a)-i matrix [Ref- 7, p. 250 ]. 



>s 



^11 ^12 
S21 S22 




(A^PA)-i or (A^A)-i = 



^ml 



’m 2 



mm 



The standard deviations of adjusted angles and 
distances can be computed by the following equation; 




1/2 



(3.25) 



where 6j is the standard deviations of the ith adjusted 
quantity. The quantity in the ith row of the 7 matrix, Q*; 
is an element of A (a"'" PA) a"'" matrix for uncorrelated 
observations of unequal precision. For uncorrelated 
observations of equal precision, Q- is an element of the 
A(a''’ A)-ia''‘ matrix [Ref. 3, p. 912]. 



A (a''’pA) -ia’’’ 
or 

A (A*^ A) -ia’’ 



^11 

Q 



21 



*n1 



■12 

•22 



•n2 



■In 



•2n 



Q 



nn 



A closed traverse was conducted at Moss Landing. 
After the second iteration, the 7 matrix can be found by 
applying the values of A, B, and X matrices to Equation 
(3.20). 



4 3 



For uncorrelated observations of unequal precision, 
tha V matrix is found to be 



^ai 


= 0.93149 


X 10-s 


radians 


or 


1.92134" 


'^a2 


= -1.63120 


X 10-5 


radians 


or 


-3.36459" 




= -1.28363 


X 10-5 


radians 


or 


-2.64768" 




= -2.30436 


X 10-5 


radians 


or 


-4.75308" 


^d1 


= 0.00024 


meters 








’d2 


= 0.00039 


meters 








^d3 


= 0.00358 


meters 









The standard deviation of unit weight is found by solving 
Equation (3.22), the result is 



(5 = ± 2.69685 



The values of PA)-^ are 



10 -s 



1.25 178 
0. 22019 
1. 25575 
0. 10788 



0.22019 1.25575 

0. 13819 ■ 0.21810 
0. 21810 1.46408 

0. 1 1127 0. 03815 



0.10788 

0.11127 

0.03815 

0.22515 



The standard deviations of the best estimate values for 
positions can be found by solving Equation (3.24), the 
results are 



^E2 


= ± 


0.00954 


meters 




= ± 


0.00317 


meters 


^E3 


= ± 


0.01032 


meters 


^M3 


= ± 


0.00405 


meters 
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The values of A ( a"’’ P A) are 



0.064 


-0. 052 


-0.021 


0.0 10 


-4.165 


-6.963 


-63.38 


-0.052 


0. 225 


-0. 175 


0.002 


3.4 15 


4.690 


41.59 


-0.021 


-0. 175 


0. 211 


-0.014 


2.563 


3. 574 


31.79 


0.010 


0. 002 


-0. 014 


0.002 


-1.813 


-1.301 


-9.50 1 


-4. 165 


3. 415 


2. 563 


-1.813 


9917 


-117.4 


-1047 


-6.963 


4. 690 


3.574 


-1.301 


-117.4 


9830 


-1525 


-63.88 


4 1.59 


31 . 79 


-9.501 


-1047 


-1525 


76323 



The standard leviations of adjusted angles and distances are 
found by solving Equation (3.25), the results are 



<531 = ± 1.40346 

dg2 = ± 2-64103 
633 = ± 2.55260 
C$34 = ± 0.26137 
= ± 0.00269 

d I 

<5^2 = ± 0.00267 
d^3 = ± 0.00745 



seconds 

seconds 

seconds 

seconds 

meters 

meters 

meters 



For uncorrelated observations of equal precision, 
the V matrix is found to be 



V 



V , = 1.26896 X 10 -s 

V 32 = -1.56893 X 10 -s 

V33 = -1.7541 4 X 10-5 

Vg4 = -2. 23358 X 10-s 

V, = 0. 14037 X 10-7 

v^^ = 0. 23844 X 10-7 

d2 

y,, = 0.2436 5 x 10-7 

d 3 



radians or 2.61742" 

radians or -3.23615" 

radians or -3.61817" 

radians or -4.60710" 

meters 

meters 

meters 



The standard deviation of unit weight is found by 
solving Equation (3.23), the result is 



=+2.01 144 X 10-5 



The values of (a'^'a)— i are 



0. 28193 0. 06827 -0.85449 0.70380 

10-» 0 . 06827 0.01664 -0.20731 0. 17075 

-0.85449 -0.20731 2.59561 -2.13724 

0.70380 0.17075 -2.13724 1.75999 

k > 

The standard deviations of the best estimate values for 

positions can be found by solving Equation (3.24), the 

results are 





0.00 107 


meters 


^N2 - - 


0.00026 


meters 


<5e3 = ± 


0.00324 


meters 


"n3=- 


0.00267 


meters 


A (A^ A) -lA^ 


are 





10-3 



1. 472 - 


27. 197 


26.960 -1.235 


-0.518 


-0. 343 


-0. 356 


-27. 197 


503 . 1 


-498. 7 22.855 


0. 348 


0. 480 


0. 473 


26.960 


-498.7 


494. 4 -22.66 


0.370 


0.523 


0.518 


-1.235 


22.855 


-22.66 1.038 


-0.20 1 


-0. 1 59 


-0. 134 


-0.518 


0. 348 


0.370 -0.201 


1 0000 


-0.001 


-0.001 


-0.843 


0.480 


0. 523 -0. 159 


-0.001 


10000 


-0.001 


-0.856 


0.473 


0. 518 -0. 134 


-0.001 


-0.001 


10000 



The standard deviations of adjusted angles and distances are 
found by solving Equation (3.25) , the results are 



^ai = ± 0.15917 seconds 

^a2 - ± 2. 94270 seconds 

^a3 = ± 2.91732 seqoads 

^^34 - ^ 0. 13370 seconds 

^di = ± 0.00002 meters 

*^d2 - * 0.00002 meters 

^d3 = ± 0.00002 meters 
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The computations in this subsection were performed 
by the NFS IBH 3033 computer using 16 decimal places and 
were rounded off to 5 decimal places prior to output. 

3 • least Squares Adjustment by, the Condition S qu at io n 
Method 



The general principle of least squares adjustment by 
the condition equation method in surveying is 



to minimiz 
squares of t 
necessary ma 
the correcti 
to zero bv a 
determined f 
from the obs 
values) ; thu 
not changed 



e a function consisting of the sum of the 
he corrections to the observations plus the 
thematical conditions involving some or all 
ons; each condition by itself is made equal 
dding the corrections to the discrepancy 
rom a preadjusted calculation (i.e. calculated 
erved values, rather than from adjusted 
s, the magnitude of the sum of the squares is 
by adding conditions [Ref. 8]. 



For uncorrelated observations of unequal precision. 



the 


pri 


ncipl€ 


! condition of 


least squares function 


may be 


exp 


ress 


ed in 


matrix 


form as 


follows : 








U = 


vTpv - 


2K[ BV + 










= 


yTpY - 


2[BV + 


W ]Tk 








= 


yTpY - 


2[ (BV)'^ 












yTry - 


2[ V''" b'*’ ♦ 


W''’ ]K 








= 


YlpY - 


2v'''b'''k 


+ 2W‘'’K 








= 


minimum 








whe 


re 0 


is th 


le least 


square 


s function matrix, K is 


a matrix 



of Lagrange multipliers, P is the weight matrix, 7 is the 
correction vector, 3 is the constant coefficients of the 
corrections, and W is the precalculated discrepancy 
[Ref. 9]. It is numerically more convenient for later 
development to multiply by 2. Taking the partial 
derivatives of the 0 matrix with respect to each of the 
corrections and equating to zero leads to 
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3 D = 2PV - 2B^K = 0 

3v 



V = P-ib''’k 



(3.2 6) 



This equation is called a correlate equation or correlate 
equation matrix- The solution of the Lagrange multipliers 
matrix can be derived by multiplying Equation (3.26) by the 
B matrix then adding the F matrix 

B7 + W = BP-ib''’k ♦ R 



(BV + W) is the condition equation matrix which must equal 
zero. The solution of the Lagrange multiplier vector is 

BP-ib''’k = -H 

K = (BP-^b"'' ) (-H) (3. 27) 

The correction vector is derived by substituting 
Equation (3.27) into Equation (3.25) 



7 = p-ib"'' (BP-ibT) - i (-») (3.2 8) 

For uncorrelated observations of equal precision, 
the correction vector can be derived similarly to the 
unequal precision case or by using the inverse of the weight 
matrix which is the identity matrix. Equation (3.28) 
becomes 



V = b’’’ (BB"’’ ) -1 (-H) 



(3.29) 



Closed traverse station positions may be adjusted by 
using the technique of least squares adjustment by the 
Condition Equation Method- There are two condition equa- 
tions. They are the azimuth and coordinate condition equa- 
tions. The coordinate condition equations are divided into 
two parts, which are E and N coordinate condition equations- 
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The azimuth condition equation is the sum of the 
corrections to the angles in a traverse plus a precalculated 
discrepancy (Equation 3.6) which must equal zero. where 
v^j equals the corrections to the angles, the general equa- 
tion can be written; 



^ V • + =0 

.si- a I 1 

i : 1 

n n -1 

2 1 t ~ d 1 



(3.30) 



where v^j is the correction to the distances. However, 
this term equals zero and does not effect the equation, but 
it does reserve space in matrix notation for distance 
corrections. A closed traverse which was conducted at Hoss 
Landing can be written with an azimuth condition of 



1. V 



at 



+ 1 . V 



a2 



1. V 



a3 



1 . V 



a4 



0. V 



dl 



+ 0. V 



d2 



0. V. 



d3 



+ Wi = 0 



(3.31) 



where the precalculated discrepancy (W^) is +8. 844" or 
+0.00004288 radians. 

The simple linearized form of an E coordinate condi- 
tion equation is the sum of the adjusted departures and must 
equal zero. 3y applying Equations (3.2) and (3-7), the 

general formula for E coordinate condition equation is 
n- 1 

2(<^i + ^ ^^af " ^ ^ (3.32) 

i r 1 

where . 

AZ^; = AZp + + V^.) - (i - 1) 1800 (3.33) 

j = 1 

AZgjis the adjusted azimuth. By using Equations (3.32) and 
(3.33) with the data conducted at Hoss Landing, these 
equations become 

(<^1 + v^,).sin AZg, + (d 2 + v^ 2 )-sin AZg 2 

- (E^ - E^ ) =0 
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<*^3 ^d3 ) ^^a3 



(3.34) 



1800 



where 



Az , 
a 1 


=r 


AZp 






^al 




= 


Az, 




^a1 






- 


AZp 






* ^a1 




=: 


Az^ 




'^al 






= 


AZp 






^ai 




= 


AZ3 




^al 


* ^a 



2 *a2 



2 *32 



3^ V33 - 3600 



AZj is the precalculated azimuth or unadjusted 
azimuth (Equation 3.1). 

From trigonometry 

sin ( A +Aa ) = sin k.cos^k + cos A.sin^A 

in which AA is a very small angle in radians, let 
cosAA = 1 and sinAA =Aa 
ther efore, 

sin { A +Aa ) = sin A +Aa.cos A 
from the first term of Equation (3.34) 

(d, + )sin AZ 21 

= (d^ + )sin (Az^ + ) 

= (d^ + v^^)[sin Az^ + cos Az^ ] 

= d^ sin Az^ + v^, d^ cos Az^ + ^ sin AZj + v^^ cos Az^ 

(3.35) 

(v^^ *^ai^ very small, by letting this term equal zero. 
Equation (3.35) becomes 

(di + v^i)sin AZ3, 

= d. sin Az, + v , d, cos Az, + v.. sin Az, (3.36) 

1 I a 1 1 1 (j 1 I 
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The second and third terms of Equation (3.34) can be derived 
similar to Equation (3.36), so 



(^2 
















= d 


2 sin AZ2 *■ 


^^al 


* 'jsjascos 


Az 


2 + 


2 


( 3 . 37 ) 


(^3 


+ v^3)sin 


■^^a 3 












= d 


3sin Az^ ■*• 


(^al ■"''a 2 -"^a 3 J ^ 3 ^°^ 


AZ3 + v^3 sin 


AZ3 


( 3 . 33 ) 




Substituting 


Equations 


( 3 . 


36 ), ( 3 . 37 ), 


and 


( 3 . 38 ) 


nto 


Equation 


( 3 . 34 ) 


and rearranging terms 






^al 


(d^ cos Az-| 


+ d2COS AZ2 + d3 


COS 


AZ3) 







a 2 


(d2COS 


AZ2 + d3COS 


Az 


3 ) 




Va3d3COS AZ3 


d1 


sin Az^ 


+ V. „sin Az 
d 2 


2 


+ 


^d3 


sin AZ3 


sin Az^ 


+ d2sin AZ2 


+ 


^3 


sin 


AZ3 — ( E j — 



This equation can be written in the general formula as 
n - 1 n -1 n -1 n -1 

2(v j * ^d.sin az. 

i rl k=i ^ i = 1 i z 1 

- (S^ - E^) =0 (3.39) 

substituting Equations (3.2), (3.3), and (3.7) into Equation 

(3.39) 



n-1 n-1 n-1 

i : 1 kz I i : 1 

For example, if the number of the observed angles (n) at 
Moss Landing is four, this equation can be expanded to 

(ANi +AN2 ^An^) + (^^2 +AN3) + 

+ sin Az^ + ^^2 * ^d3^^^ ^^3 * ^2 ~ 

5 1 



0 



Substituting the numerical values from the precalculated of 
this traverse into this equation, the result is 

(6702.37612) + (53 18.48345) V 32 + (4 999 . 28 284) v^^ + °*^a 4 

-(0.23567) +{0.48818)7^2 + (0 - 63 570) ^ + »2 ^ ° (3-41) 

where the precalculated discrepancy (tr^) is 0-08608 m and 
O.v^^is equal to zero, which does not effect this equation, 
but space is reserved in matrix notation for the last angle 
correction. 

The N coordinate condition equation is the sum of 
the adjusted latitudes and must equal zero. By applying 

Equations (3.3) and (3.7), the general condition equation is 
n-1 

■ (N-r ■ N,) =0 

I X 1 

This equation can be derived similarly to the E coordinate 

condition equation case, resulting in 
n - 1 n - 1 n-1 

Az. + W 3 = 0 (3.4 2) 

iri k^l i=l 

From the data acquired at Moss Landing, this equation can 
be expanded to 

-v^^ (Ae, +AE 2 ^Ae^) - V 32 (AE 2 ^Ae^) + V 33 AE 3 

+ cos AZ- + V cos AZ -5 + V. cos Az, + H -1 = 0 

d 1 1 d 2 2 d3 3-5 

Substituting the numerical values from the calculation of 

this traverse into this equation, the result is 

-(3959.89460) v_ - (4295. 49626) V- 7 ( 4 1 1 6 . 94 755) v _ + O.v^ 

a I 3 c 3 ^ 3^ 

+ (0. 97183) V + (0. 87274) v^^ + (0 . 77 1 94) v = 0 (3.43) 

(j 1 d 2 Q J J 

where the precalculated discrepancy (W^) is -0.08936 m. 
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Equations (3.30), (3.40), and (3.42) can be written 

in the matrix form as (BV + W) , where 3 is the constant 
coefficients of the corrections, 7 is the correction vector 
which is composed of the angle and distance corrections, W 
is the precalculated discrepancy. Equations (3.31), (3.4 1), 

and (3.43) are written in matrix form as 



✓ 










1 


6702.37612 


-3959.89460 




’al 


1 


5318. 48345 


-4295.49626 




^a2 


1 


4999. 28284 


-4116.94755 




^a3 


1 


0.00000 


0. 00000 


7 = 


^a4 


0 


-0.23567 


0.97183 




^di 


0 


0.48818 


0. 87274 




^d2 


0 


0.63570 


0. 77194 




^d3 

< > 






✓ 





-H 



-W, = -0. 00004 

= -0- 08608 
-ff- = 0.08936 

N J J 



For uncorrelated observations of unequal precision, 
the inverse of the weight matrix is the diagonal matrix, 
which is equal to the observation’s variance. The diagonal 
elements of the inverse of the weight matrix is 



p -1 

11 


— 


sin2 (1. 984”) 


— 


0.92519 


X 


10-10 


p-i 

2 2 


= 


sin2 (1.405”) 


= 


0. 46398 


X 


10-10 


p -1 

33 


= 


sin2 (1. 203”) 


- 


0.34016 


X 


10-10 


p-l 
4 4 


= 


sin2 (1.614”) 


- 


0.61229 


X 


10-10 


p-l 

5 5 


= 


(0.001) 2 


= 


0.10000 


X 


10-5 


p-l 

66 


= 


(0.001) 2 




0.10000 


X 


10-5 


p-l 

7 7 


= 


(0.003) 2 




0.90000 


X 


10-5 
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For uncorrelated observations of unequal precision, 
the correction vector is found by solving Equation (3,28) 





’a1 


= 0.93141 


X 10-5 


radians 


or 


1.92117” 




CM 


= -1.63114 


z 10-5 


radians 


or 


-3. 36446” 






= -1.28358 


X 10-5 


radians 


or 


-2.64758” 


V = 


’a 4 


= -2.30433 


X 10-5 


radians 


or 


-4.75312” 




^d1 


= 0.00024 


meters 










^d2 


= 0.00039 


meters 










a 


= 0.00357 


meters 









The corrected angles are 



ai 


= 2460 


05 ' 


45. 12117” 


a2 


= 2220 


51 • 


05. 23554" 


CC3 


= 1900 


14 • 


59. 95242" 


^4 


= 2770 


05' 


12. 24688" 



The corrected distances are 

di = 1424. 00424 m. 
d.2 = 365. 74439 m. 

d 3 = 6476.27457 b. 

Using the corrected data to recompute the coordinates is 
similar to the initial traverse computations. The adjusted 
coordinates are 

E 2 = 607,943.45521 m. 

H 2 = 4,073,939.74809 m. 

E 3 = 608,122.00139 m. 

N 3 = 4,074,258.95029 m. 

E^ = 612,238.85256 m. 

= 4,079,258.31754 m. 

The last adjusted coordinates are the same values as the 
fixed coordinates. Hence, it is not necessary to iterate. 
For uncorrelated observations with equal precision, the 
correction vector is found by solving Equation (3.29) 
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V = 



^al 


= 1.26884 


X 


10-5 


radians 


or 


2. 61 717" 


^a2 


= -1.56886 


X 


10-5 


radians 


or 


-3. 23601" 


V 


= -1.75408 


X 


10-5 


radians 


or 


-3. 61804" 


a3 














^a4 


= -2. 23360 


X 


10-5 


radians 


or 


-4.60713" 


^di 


= 0. 14034 


X 


10-7 


meters 






C\J 


= 0.23842 


X 


10-7 


meters 






^d3 


= 0.24362 


X 


10-7 


meters 







The corrected angles are 





= 


2460 05’ 45.81717" 






= 


2220 51* 05.36399" 




C<3 


= 


1900 14' 58.98196" 






= 


2770 05' 12.39287" 




The corrected distan 


ces 


are 






= 


1424. 00400 m. 




CM 


= 


365.74400 m. 






= 


6476.27100 m. 




By using the corrected 


data to recompute 


the 


the adjusted coordinates are 




^2 


= 


607,943.45994 


m • 




= 


4,073,939.74899 


in. 


^3 




608, 122.00770 


m • 






: 4,074,258.95013 


m. 


^4 




612,238.85256 


m. 


»^4 




4,079,258.31754 


ni. 



The last adjusted coordinates are the same values as the 
fixed coordinates. Hence, it is not necessary to iterate. 
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The standard deviation of an observation which has 
unit weight can be found by the following equations. For 
uncorrelated observations of unequal precision. 



(5- = 



v^pv 



1/2 



(3.44) 



and, for uncorrelated observations of equal precision. 



<5 



o 






1 /2 



(3.45) 



r 

where r is the number of condition equations. There are 
three condition equations for a closed traverse (r=3) . The 
standard deviations of adjusted angles and distances can be 
found by the following equation: 





1 /2 



(3.46) 



where (5j is the standard deviation of the ith adjusted 
quantity. The quantity in the ith row of the 7 matrix, Q|j 
is an element of the P-^ - P b'*' (BP b"'") BP-^ matrix for 
uncorrelated observations of unequal precision 
[Ref. 3, p. 918]. For uncorrelated observations of equal 
precision, Qjj is the element of the I - B^ (Bb"^) -^B matrix. 
The matrix P-i - p-i b"^ (BP-^ B*^ ) — ^B?-i or the matrix 
I - b"*" (Be"^) — IB is equal to 



11 


Q 1 2 




• 


• 


• 


• 


Qln 


21 


Q 22 




• 


• 


• 


• 


220 



-n1 






• 2 nn 
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For uncorrelated observations of unequal precision, 
the standard deviation of unit weight is found by solving 
Equation (3.44), the result is 



6 = ± 2.69679 

o 



The results of P- 


-1 - p - 


13 ”^ (BP-ib''’) -1 3? 


-1 are 








0.064 


- 0. 052 


- 0.021 


0.0 10 


- 4.165 


- 6.963 


- 63.68 




- 0. 052 


0 , 225 


- 0. 175 


0.002 


3.4 15 


4.689 


41 .59 




- 0.021 


- 0. 175 


0. 211 - 0.0 14 


2 . 563 


3.574 


31.79 


10-10 


0 . 010 


0 . 002 


- 0.014 


0.002 


- 1.813 


- 1. 301 


- 9. 500 




- 4.165 


3-415 
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The standard deviations of adjusted angles and distances are 
found by solving Equation (3.46) ; the results are 



«al 




± 


1.40340 


seconds 




= 


± 


2. 64096 


seconds 


<5.3 


= 


± 


2.55253 


seconds 


<5.4 


= 


± 


0.26136 


seconds 


^dl 


= 


± 


0.00269 


meters 


^d2 


= 


± 


0.00267 


meters 


^d3 


= 


± 


0.00745 


meters 



For uncorrelated observations of equal precision, the 
standard deviation of unit weight is found by solving 
Equation (3.45), the result is 



6 , = ± 2-01139 X 10-5 
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The results of I 



b''’ (BB’’’ )— 1 B 



are 
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The standard deviations of adjusted angles and distances are 
found by solving Equation (3.46) ; the results are 



<^al 
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± 


0. 15918 


seconds 


^a2 
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± 


2.94262 


seconds 


-3a 3 


= 


± 


2.91722 


seconds 


-3a 4 


= 


± 


0. 13369 


seconds 




= 


± 


0.00002 


meters 




= 


± 


0.00002 


meters 


3.3 


= 


± 


0.00002 


meters 



The computations in this subsection were also 
performed on the NPS IBM 3033 computer; they used 16 decimal 
places, and were rounded off to 5 decimal places in the 
final solutions. These are more decimal places than are 
normally used in practice because such precisions are not 
generally attainable by corresponding observations. 

However, in comparing two computational methods that are 
theoretically equal, one method may be more sensitive to 
round off error than the other. This will be commented on 
in the next chapter. 
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IV. DISCUSSIONS AND AMiilSIS OF RESULTS 



Three programs were used to compute the initial traverse 
and adjust the closed traverse station positions. Programs 
1, 2 , and 3 were used to compute and adjust the traverse 
station positions by the Approximate Method, the Indirect 
Observations iMethod, and the Condition Equation Method, 
respectively. The programs were written in KATFIV language 
for implementation on the NPS IBM 3033 computer. The 
computer output and listings of Programs 1, 2, and 3 are 
provided in Appendices A, B, and C, respectively. The 
maximum number of intermediate traverse station positions 
that can be computed and adjusted by these programs is 30. 
The programs were tested using several fictitious data sets 
to ensure their performance in handling the various 
intermediate traverse station positions. 

The computer storage area and CPU time of Programs 1, 2, 
and 3 has been compared (Table X) . The adjustment of a 
closed traverse by the Approximate Method did not use a 
weight matrix in the matrix computations. The computer 
program for this method was written by using the variables 
in only one dimension for storage of both data and results. 
This method did not require any iterations. The least 
squares adjustment of a closed traverse characteristically 
is used to simultaneously eliminate closing errors in 
azimuths and distances, and Programs 2 and 3 both utilized 
such an adjustment. The computation by either Least Squares 
Method used matrices in two dimensions to compute the 
correction vector. Likewise the computer programs were 
written using variables in two dimensions. Included are all 
subroutines from Program 1 and extra subroutines for each 
individual program. Therefore, both Programs 2 and 3 used 
more computer storage area and CPU time than Program 1. 
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TABLE X 

The Comparisons Between the Computer Storage Area and CPU 

Programs 2 , and 3 



M-l 
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0 > 

B 

•H 

EH 
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The Indirect Observations ilethod requires that the 
number of observation equations be equal to the sum of 
observation equations of both observed angles and distances. 
The number of observation equations is not constant. It 
changes depending on the number of the traverse stations-- 
this means that the row dimension of each matrix will change 
depending on the number of observation equations. 

Therefore, the subroutines affected by the number of trav- 
erse stations are difficult to write in the general program. 
The coefficients in the A matrix are computed by taking- the 
partial derivatives of the functions with respect to each 
unknown variable, where the number of unknown variables is 
not constant. The number of unknown variables changes 
depending on the number of the traverse stations. 

Therefore, the subroutine for the computed A matrix is also 
difficult to write in the general program. The Condition 
Equations Method has only three equations which makes it 
easier to derive the general form and write the computer 
program. The adjustment using the Indirect Observations 
Method does not apply the corrections to the observed values 
directly; therefore, it requires the approximate values of 
the unknown coordinates be computed for the correction 
vector of the unknown coordinates. This method requires at 
least two iterations to check the insignificance of the 
correction vector when it is compared with an arbitrarily 
selected small number. The adjusted coordinates from the 
Condition Equations Method can be checked at the first iter- 
ation. Thus, Program 3 is more economical as it uses less 
computer storage area and CPU time than Program 2 (Table X) . 

The closed traverse at Moss Landing originated and 
terminated at control points with known positions which were 
determined by third-order methods. The azimuth closure was 
2.211" per station and the position closure was 1:66,617, 

The classification, standards of accuracy, and general 
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specifications for a third-order class I traverse (Table I) 
indicate the azimuth closure is not to exceed 3" per station 
and position closure must be better than 1:10,000. This 
traverse met the specifications and standards of accuracy 
for a third-order class I traverse. 

Using Programs 1,2, and 3, the traverse was computed 
and adjusted in UTH grid coordinates. The final adjusted 
coordinates were transformed from OTM grid coordinates to 
geographic coordinates [Ref. 4, pp. 319-321]. The differ- 
ence between the coordinates for each method and NOS 
[Ref. 10] have been computed (Table XI). The technique of 
Least Squares for both methods yields identical computa- 
tional results to at least five decimal places (Table XI, 
Methods 3.1 and 4.1; 3.2 and 4.2). Least Squares provides 

the best estimates for positions of all traverse stations. 
Program 2 performs a statistical test, yielding the standard 
deviations of both adjusted positions and observed values. 
Conversely, Program 3 only yields the standard deviation of 
adjusted observed values.- Comparisons of the standard devi- 
ations of adjusted observed angles and distances were made 
between these computed by the Indirect Observations Method 
and the Condition Equations Method (Table XII) . The 
computed standard deviations of the adjusted angles differed 
significantly in the fourth decimal place in several cases, 
which is a slight indication that the Indirect Observations 
Method is more sensitive to round off error. The standard 
deviation estimates are larger for the Indirect Observations 
Method in every case of a significant difference in the 
fourth decimal place. However, a significant difference in 
the fourth decimal place is insignificant in terms of the 
observational precision. The standard deviations of the 
adjusted distances for both methods yield identical results 
(Table XII). Both Methods differ from the NOS positions in 
the fourth decimal place in seconds of arc for both latitude 
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and longitude, 
dards recuires 
of arc for both 
can be used for 



The horizontal control for third-order stan- 
accuracy to three decimal places in seconds 
latitude and longitude--thus, these methods 
computing third-order positions. 
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Geographic Coordinates 
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Unequal precision 36 48 35.38122 -0.00009 121 47 16.39152 +0.00015 



The Comparisons Between the Standard Deviations of Least Squares 
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7, CONCIDSIOMS AND RECOMMENDATION 



A. CONCLUSIONS 

The conputei: programs developed by the author and 
contained in this -thesis have a variety of useful and prac- 
tical hydr o<sra phic applications. In hydrography, geodetic 
field work plays a key role. Horizontal control is neces- 
sary in determining positions and, even in areas which 
appear to initially have adequate control, additional 
control stations often need to be established after arrival 
in the field. These field-generated control points must ,:e 
adjusted to fit '.-'ithin the existing survey net and to ain:.- 
mize errors in the field neasur-oaents. The errors can he 
minimized by these computer programs. Direct application of 
these programs would be enormously beneficial, since no 
equivalent software exist at NPS to perform such adjustments 
at the present time. 

In the Approximate Method, a traverse is adjusted by 
computations using a hand calculator. This method is suit- 
able for field computation and the adjusted coordinates meet 
the specifications and standards of accuracy for a third- 
order class I traverse. Although the accuracy from this 
method is less than the Least Squares Metnod, it does not 
require a computer and it furnishes coordinates which can be 
checked in the field by the field party. 

The Least Squares Method provides the best estimates for 
positions of all traverse stations. Least Squares does 
require more of a computational effort than the Approximate 
Method. However, the accuracy of Least Squares is better 
than the Approximate Method, and the Least Squares Method 
should be used for the final office computations. Both 
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Least Squares Methods yield identical results, but the 
Condition Equations Method is more economical and easier to 
derive than the Indirect Observations Method. Therefore, 
the Condition Equations Method should be applied at MPS. 

B. EECOSHENDATIOH 



Studies should be continued at NPS to compare the 
economy of adjustment methods used in this thesis with other 
methods of least squares adjustment. 
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